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I.  INTRODUCTION 


The  compressible  flow  boundary'  layer  program,  COMPBL^ ,  was  written  to 
calculate  the  two  dimensional  boundary  layer  flow  for  an  arbitrary  pressure 
gradient  in  air.  An  integral  technique  based  on  the  WalK^  method  was  utilized 
in  this  program.  Capabilities  include  axisymmetri c/plane  geometry, 
laminar/turbulent  flow,  incompressible/compressible  flow,  with  or  without 
heat  transfer.  This  program  was  written  for  the  CDC  CYBER  class  of  computers 
and  is  in  use  locally  in  the  BATCH  mode. 

There  are  two  deficiencies  of  this  code  as  it  presently  is  used: 

A.  It  is  limited  to  boundary  layer  flows  in  air. 

B.  It  has  limited  usefulness  as  an  aero  design  tool  since  it  is  uti¬ 
lized  in  a  BATCH  mode. 

These  two  deficiencies  w^re  correcied  in  this  analysis  and  insulted  in 
COMPBL-II  which  has  the  following  additional  capabilities: 

A.  Allows  an  arbitrary  chemical  composition  of  the  boundary  layer 

flow. 

B.  Runs  either  interactively  on  the  INTERDATA  8/32  or  BATCH  on  the 
CYBER  Series. 

II,  ANALYSIS 

The  chemical  species  chosen  for  use  in  this  analysis  were  taken  from  the  DATABANK 
subroutine  of  the  JANNAF  Standardized  Plume  Plowfield  (SPF)  Code.  These  include 
those  of  most  interest  for  rocket  calci’-Iatiori.  It  was  necessary  to  utilize 
transport  property  data  for  the  chemical  species  chosen  in  orde^  to  calculate 
the  boundary  layer  for  internal  flows.  These  data  were  tabulated  over  a  wide 
range  of  temperatures  (200-3700  °K)  covering  both  internal  and  external  flows. 
These  data  include: 

A.  Elemental  viscosity,  Hi 

B.  Thermal  conductivity,  ki 


C.  Constant  pressure  specific  heat,  Gp^ 


Data  a  and  b  were  taken  directly  from  NASA  TR-R-132^.  Data  c  was  taken  from  the 
DATABANK  subroutine  of  the  SPF  code  and  include  coefficients  of  a  curve  fit  of 
Cp^  according  to  the  relation 


Q 

Pi  =  A  +  BT  +  CT^  +  DT5  +  ET4 
R 


(1) 
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These  data  were  assimilated  into  mixture  properties  assuming  ideal  gas  beha¬ 
vior  for  both  thermal  conductivity  and  specific  heats. 

Ic  .  =  2  vl  k. 


•^mix  i=l  i 


=  I 


For  the  mixture  visccsity,  this  simple  summation  is  inadequate  and  a  method 
for  the  determination  of  mixture  -viscosities  found  in  Walker^  was  utilised. 
Because  of  the  property  changes  from  air  to  an  arbitrary  chemical  mi::ture,  it 
was  necessary  to  account  for  changes  in  the  laminar  recovery  factor,  turbulent 
recovery  factor,  laminar  Reynolds  analogy  factor,  and  the  turbulent  Reynolds 
analogy  factor. 

The  -viscosity  power  law  assumption  was  utilized  in  tnis  investigation  -*-0 
determine  the  exponent,  w,  for  use  in  the  boundary  layer  calculations.  This 
assumption  is  valid  when  £ny  is  linear  with  UnT  as  shown  in  Figure  1 . 


Figure  1.  Viscosity  as  a  function  of  temperature. 

This  linearity  assumption  was  made  in  the  program  and  the  method  of  least 
squares  was  utilized  to  determine  the  slope  w. 

u  /  T  \ 

ho  '^o  /T\ 

£np  =  JlnVo  +  win  1-:;^) 

TV  '  ^O' 

^  ^wJn  +  c 


(-t) 


Equation  (5)  was  curve  fit  using  a  linear  least  squares  routine, 
ficients  c  and  w  were  evaluated  utilizing  (6)  and  (7). 


The  coef- 


ax^  -  noxx 
Ov  -  box 


The  derivation  of  (6)  and  (7)  is  shown  in  Appendix  A. 


III.  THE  CODE 


Subroutine  NOAIR  is  the  principal  addition  to  COMPEL.  If  the  program  is  run 
for  a  gas  mixtare  other  than  air,  Subroutine  NOAIR  is  utilized. 

Inputs  into  NOAIR  include  the  Mach  number  and  pressure,  the  elements 
comprising  the  mixture  (along  with  the  mole  fraction  of  each),  and  the  mixture 
temperature. 

From  these,  the  mixture  viscosity  is  calculated  at  the  input  temperature  by 
using  the  cubic  spline  interpolation  routine  mentioned  before.  Using  this 
■viscosity,  the  viscosity  power  law  exponent  is  derived  by  using  the  least 
squares  fit  routine.  Also,  ;he  specific  beat  ratio  is  found  using  the  equation; 


^Pmlx  “  ^ 


The  Prandtl  number  is  also  returned  by  the  equation: 


''  -{-f-)  (1^)  (tBc) 


The  Prandtl  number  is  then  used  to  change  the  laminar  and  turbulent  recovery 
factors,  as  well  as  the  laminar  and  turbulent  Reynolds  analogy  factors: 


RT  =  Pr 


SL  =  Pr 
ST  =  Pr 


The  density  and  velocity  of  flow  also  are  calculated  as  -Tollows; 


pMWmix 


[s.ElAS  ] 


u  —  M  V 


These  are  used  to  calculate  the  Reynolds  number. 


„  pu 
Re  =  — 
P 


A.  Data  Tapes 

This  program  is  designed  to  run  on  the  CYBER  74  as  well  as  run  interac¬ 
tively  on  the  INTERDATA  8/32.  It  is  written  in  Standard  Fortran  IV. 

Both  versions  contain  three  (3)  DATA  FILES  within  Subroutine  NOAIR; 

C1:  Contains  individual  species  viscosities  for  temperature  range 
200-3700*K. 

COEF;  Contains  constant  pressure  specific  heat  coefficients  taken  from 
Program  SPF  DATABANK.  (Used  lo  calculate  Cp^/R) 

XKAY:  Contains  individual  species  thermal  conductivity  for  temperature 
range  200-3700 “K. 


B.  Arrangement  of  Data 

To  make  Subroutine  NOAIR  function  as  simply  as  possible,  it  was 
necessary  to  arrange  the  data  in  a  logical  order.  From  the  listing  of  elements, 
it  is  evident  that  the  elements  are  in  semi-alphabetical  order,  and  are  numbered 
as  such. 

These  numbers  are  important,  because  the  arrangement  of  the  data  (WTMOLE, 

Cl,  COEF,  XKAY)  also  follow  this  order.  In  other  words,  for  element  N2,  the 
corresponding  number  is  31*  Now,  the  molecular  weight  for  N2  is  the  31st  ele¬ 
ment  in  the  array  WTMOLE.  The  other  three  data  files  are  also  arranged  in  this 
way;  however,  in  these  latter  data  files,  there  are  four  lines  of  data  per  ele¬ 
ment  for  Cl  and  XKAY,  and  two  lines  of  data  per  element  in  COEF,  e.g.,  take  N2, 
number  31  •  The  viscosities  for  temperature  range  200-3700*K  would  be  lines 
121-124  in  arrays  Cl  and  XKAY  and  lines  61-62  in  COEF  (See  Table  I). 

C.  Adding  New  Elements 

Taking  the  previous  section  into  account,  it  is  necessary  to  maintain 
this  order.  To  add  new  elements,  it  should  be  easier  to  add  the  new  element  in 
position  47,  and  the  data  for  WTMOLE,  Cl,  COEF,  XKAY  should  be  attached  to  the 
end  of  each  file. 

Adding  the  new  element  will  be  c.  little  different  on  the  8/32  and  the  CYBER 
74.  On  the  8/32,  the  element  will  be  aided  at  the  end  of  the  listing: 


On  the  CYBER  74 »  the  new  element  will  be  added  on  the  end  of  the  DATA  Statement 
for  ILEM: 

DATA  (ILEM(I),  I  =  1,50)  2HAK,  3HBCL,  4HBCL2,  4HBCL3, 


$1H0,  2H02,  2H03,  2HOH,  2HNE/ 


and  the  number  47  will  be  added  to  KLEM; 

DATA  (KLEM(K),  K=1 ,50)/1 ,2,3,4,5,6,7,8,9, 


$44,45,46,47,3*0/ 

IV.  DATA  INPUT  -  CYBER  BATCH  VERSION 

CARD  1 ;  This  card  tells  whether  the  boundary  layer  program  is  to  run  on 
properties  of  air  or  of  another  gas  mixture 

In  Column  5,  type 

1  -  for  air 

2  -  for  another  mixture 

CARD  2:  (Only  if  '2*  is  typed  on  CARD  l) 

Right  Justify; 

Columns  1-10;  Mp''h  or  M*  number 
Columns  11-20;  Pressure  (in  kilopascals) 

(NOTE;  Make  sure  pressure  is  entered  in  kilopascals) . 

CARD  3:  (Only  if  '2'  is  typed  on  CARD  l) 

Number  of  elements  in  the  mixture 
Right  justify  in  columns  1-5. 
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CARD  4t  (Only  if  '2'  ia  typed  on  CARD  1) 

These  next  cards  contain  information  about  the  mixture. 

Left  justify  in  columns  1-6 

Type  in  the  element  name 

Right  justify  in  columns  10-20 

Type  in  the  mole  fraction  of  that  element 

(NOTE:  Every  element  -  mole  fraction  association  must  be  on  a  separate  card. 

There  should  be  as  many  of  these  cards  so  to  equal  the  number  typed  on 
Card  3.  Also,  mole  fractions  must  add  up  to  equal  1.0.) 

CARD  5:  (Only  if  *2*  is  typed  on  CARD  l) 

Right  justify  in  columns  1-10 
Temperature  of  the  mixture 

(NOTE:  Temperature  must  be  entered  as  degrees  Kelvin  and  the  acceptable  range 
is  200-3700 ‘K.) 

CARD  6:  Right  justify  in  Columns  1-80 

Place  the  title  of  the  boundary  layer  run.  It  is  advised  that  this 
title  be  centered. 

CARD  7j  This  card  contains  information  necessary  for  boundary  layer 
computations; 

All  are  left  justified 

Columns  1-10:  Enter  the  flow  type  (Laminar  or  Turbulent) 

Columns  11 -20:  Enter  geometry  type  (Plane  2-D  or  Axisym) 

Columns  21 -30:  Enter  velocity  input  data  type  (Mach  or  Mstar) 

Columns  31-40:  Enter  wall  temperature  data  (Theta  or  Twtinf) 

CARD  8:  This  card  is  used  to  override  the  default  values  in 

Namelist  BL. 

Place  a  $  in  column  2  and  the  name  of  the  namelist,  type  the  name  of 
the  variable  along  with  its  new  value,  separate  all  variables  by 
commas,  end  with  $. 
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Example: 


1'" 

2 

$ 

BL  HSTART 

1.6,  MINF  =  0.308,  G  =  1.4$ 


CARD  9'  The  rest  of  the  cards  are  data  points  necessary  for  boundary  layer 
computations,  as  described  in  Dutton  and  AddyV ' K 

V.  DATA  INPUT  -  INTERDATA  8/32 


In  contrast  to  the  CYBER  Batch  Version  the  INTERDATA  8/32  Version  is  interac¬ 
tive.  The  user  utilizes  a  prompting  system  to  input  the  data.  The  list  of 
prompts  and  a  description  of  each  input  follows: 

1 .  Is  this  being  run  on  a  Tektronix  Terminal? 

Enter  .  Y  .  OR  .  N  . 

This  question  has  to  do  with  the  terminal  the  program  is  run  on.  If  "Y" 
is  typed  in,  the  screen  will  automatically  clear  itself  when  it  is  full. 

2.  Is  this  run  to  be  based  on  the  qualities  of  air?  Enter  Y  or  N. 

If  "Y"  is  entered,  then  Subroutine  NOAIR  is  by- passed  and  the  program 

utilizes  the  properties  of  air. 

(NOTE:  Questions  3-7  will  appear  ONLY  if  *’N*'  is  entered  for  question  ) 

3*  Please  enter  Mach  or  Mstar  Number. 

Enter  the  Mach  or  M  number.  Be  sure  to  include  decimal  point 
Ex:  >  0.9667 

4.  Please  enter  pressure  in  kilopascals 

Enter  the  pressure  in  kilopascals.  Be  sure  to  include  decimal  povnt. 

Ex:  >  20.  or  20.0 

5.  Now  enter  the  number  of  elements  in  the  gas  mixture.  The  number  can  be  bet- 

tween  one  and  99.  If  the  number  entered  is  one  digit,  include  leaning 

zero. 

Ex:  >  03  ;  >  10 
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6.  Here  Is  a  list  of  elements: 


1-AR 

2-BCL 

3-BCL2 

4-BCL3 

5-BF 

6-BP2 

7-BF3 

8-BO 

9-3203 

10-C 

11 -CH 

12-C2H2 

13-CH4 

14-C2H4 

15-CO 

16-C02 

17-CL 

18-CL2 

19-CLF 

20-CLF3 

21 -F 

22-F2 

23-H 

24-H2 

25-H20 

26-H202 

27-HCL 

28-HP 

29-HE 

30-N 

31-N2 

32-NF3 

33-NH 

34-NH3 

35-NO 

36-NO2 

37-NOCL 

38-N20 

39-NA 

40-NACL 

42-NAO 

42-NAOH 

43-0 

44-02 

45-03 

46-OH 

Enter  the  corresponding  nuir.ber(s)  of  the  elements  to  be  used  along  with 
the  mole  fraction  of  each  (i.e.,  35$  of  03»  enter  45i0.35) 

These  are  all  the  elements  contained  in  Subroutine  NOAIR.  Enter  the  element 
number,  not  the  element  name,  along  with  the  mole  fraction  of  that  element. 

Each  prompt  should  contain  the  element  number,  followed  by  a  comma  and  the  ele¬ 
mental  mole  fraction.  There  will  be  as  many  prompts  as  the  number  entered  for 
question  5» 

Ex:  for  3556  of  O3:  >  45,  0.35 

(Mole  fractions  must  add  up  to  equal  1.0.) 

7.  Now  enter  the  temperature  of  the  mixture. 

Enter  the  temperature  of  the  mixture  in  degrees  Kelvin.  Be  sure  to  include 
decimal  point. 

Ex:  >  300.00 

(Temperature  must  be  within  the  range  of  200-3700'Kelvin. ) 

The  following  cues  appear  regardless  of  the  answer  to  question  2; 

8.  Enter  type  of  flow; 

Enter. .. .LAMI. .. .for  Laminar 
Enter. . . .TURB. . . . for  Turbulent 

Variable  that  describes  the  initial  flow  regime  of  boundary  layer.  Type  in  LAMI 
for  laminar  flow  or  TURB  for  turbulent  flow. 
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9.  Sntfer  type  of  geometry: 

Inter. .. .2-D. ... fcr  Plane  2-D 
■■''nter. . .  .AXl.S. . .  .for  , Vi  symmetric 

Variable  that  describes  the  geometry.  Enter  2-D  for  plane  2-D  or  AXIS  for 
axisymmetric. 

10.  E.itev  velocity  data  input  type: 

Enter. .. .MACH. .. .for  Mach 
Enter.... MSTA.... for  Mstar  (H^) 

Variable  that  describes  whether  velocity  input  data  is  in  terms  of  Mach  number 
or  M*  (=u/c*).  Enter  MACH  for  Mach  number  or  MSTA  for  M*. 

11.  Enter  wall  temperature  input  data: 

Enter. .. .THET. .. .for  Theta 
Enter. ...TWTI.... for  Twtinf 

Variable  that  describes  whether  wall  temperature  input  data  is  in  terms  of 
d  =  ( (Te-Tw)/Te-TA)J  or  Tw/T«  .  Enter  THET  for  theta  (e)  or  TWTI  for 
Tw/Tw  .  (There  are  two  prompts  for  the  next  question.  Only  one  will 
appear  and  it  depends  on  the  answer  to  question  2). 

12.  Tne  values  for  Mach  or  Mstar  number,  Reynolds  number,  Gamma, W,RL,RT,SL, ST,  have 
beec  calculated.  Here  are  the  rest  of  the  input  variables: 

ZSTART  =  0.0000 

H3TART  =  1.572 

3START  =  0.0000 

XTRAN  =  0.0000 

?ST  =  0.0000 

EPS  =  0.1000S-05 

ERROR  =  F 

A'OTRAN  =  F 

HTCORR  =  F 
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Do  you  i/lah  to  ch&nge  any  of  the  above? 

Enter. . . .Y. . . .or. . . .N. *. . 

This  version  appears  if  the  sncwer  to  question  2  is  "N". 

Here  are  the  input  variables  with  their  present  (default)  values 

ZSTART  =  0.0000 

HSTART  =  1.572 

BSTART  =  0.0000 

MI  =  0.0000  (MINF) 

REINFL  =  0.0000 

XTRAN  =  O.OCOO 

FST  =  0.0000  (FSTINT) 

X  =•  1.405  (G) 

Vf  *  0.7000 

RL  »  0.8500 

RTR  =  0.8800  (RT) 

SL  =  0.8000 

STS  =  0.8200  (ST) 

EPS  =  0.1000E-05 

ERROR  =  F 

NOTRAN  =  F 

HTCORR  =  F 

Do  you  wish  to  change  any  of  the  above? 

Enter. . . .Y. . . .or. . . .N. . . . 

This  version  appears  if  the  answer  to  question  2  is  "Y". 
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In  either  case,  the  question  Is  designed  to  enable  the  user  to  override  any  o: 
the  default  values  for  these  variables.  (See  Figure  2  for  explanation  of  the 
variables.)  If  new  values  are  to  be  entered,  type  "Y”;  if  not,  type  "N". 

Question  13  will  appear  only  if  answer  to  12  is  "Y". 

(Two  variables  must  be  changed  from  0.0  to  another  number;  i.e.,  they  must  be 
given  a  value.  These  are  MI  (MINF)  and  REINFL,  and  need  be  changed  only  if 
question  12B  appears.) 

13.  Entfrr  names  and  new  values  (Name  =  New  Value) 

Enter  them  one  at  a  time,  pushing  return  after  each  or  enter  .  ^  .  to  signal  e 
more  variables  to  be  changed. 

> 

>% 

Enter  the  variable  name  followed  by  an  equals  sign  and  the  new  value.  To  ter¬ 
minate  changes  to  these  variables,  type  in 

Ex;  >  HSTART  =  1.6 

>  MI  =  1.0 

>  % 

VI.  SAMPLE  INPUT 

Sample  Card  Deck  Air  Case 


1 

Tdhether  run  on  air  or  not 

1  -  Air 

GS 

R3  BOUNDAI 

tY  UY5P 

Title  of  run 

TURBULENT 

AXISYM 

MACH 

THETA 

Flow,  geometry,  velocity  input  data,  wall 

temperature  data 

$BL  HSTART= 

=1  .6,MINF=0.9667,REINFL=0. 

1081E6, 

Namelist  BL 

G=1.39'i 

0.0 

1.888 

1.0 

1.888 

Local  data  cards 

1.027 

2.256 

1.707 

2.256 

2.197 

2.792 

2.11 

2.792 

5.09 

3.145 

2.306 

3.145 

XI,  YI,  MDI,  RI 

4.3 

3-5465 

2.488 

3.546‘. 

5.222 

3.802 

2.589 

3.802 

THETI 

6.157 

4.012 

2.666 

4.012 

6.63 

4.1 

2.696 

4.1 

1  5  11  21  31  41  51 


13 


ZSTAHT--SrARnNG  VALUF  OF  Z,  fLJ  f DFPAULT*0,0) 

HSTART— STA«TI-J6  VALUF  OF  H»  SKF  WALZ  PP  267-268  AND  FIGS  I.l 
AND  6.5  AND  ACCOMPANYING  REPORT.  TABLES  I  AND  II 
AND  FIG.  2  (DEFAULTS!. 57?) 

BSTART— aSTART^PI  IS  THE  INCLUDED  ANGLE  AT  THE  LEADING  EDGE 
FOR  BOTH  PLANE  2-D  AND  AKISYMMETRIC  EXTERNAL  FLOWS 
(DEFAULTaO.O) 

MIMF— — MACH  NUMBER  (OR  HSTAR)  OF  APPROACH  FLOW 
8KINFL—8EYN0L0S  NUMBER  DIVIDED  BY  CHARACTERISTIC  LENGTH  OF 
approach  flow,  l.E.  RH0INF*UTMF/MUINF.  ILJ**(-1) 
ITRAN— X  LOCATION  FOR  SPECIFIED  TRANSITION,  (LJ 
(OEFAULTsO.O) 

rSTiNT— FREESTREAM  TURBULENCE  INTENSITY— USED  IN  TRANSITION 
SUBROUTINE,  IN  PERCENT  (DEFAULTsO.O) 

RATIO  OF  SPECIFIC  HEATS  (DEFAULTS!  .405) 

H.....— exponent  on  viscosity  power  law  (DEFAULT30.7) 
HI.-— laminar  RECOVERY  FACTOR  ( DEFAULTsO. 85) 

BT««— •turbulent  recovery  FACTOR  (DEFAULTsO. 88) 

•L—— laminar  REYNOLDS  ANALOG  FACTOR  (DEFAULTsO, 80) 

ST— •••TURBULENT  REYNOLDS  ANALOG  FACTOR  (DEFAULTsO. 82) 
.CONVERGENCE  CRITERION  VARIABLE  (DEFAULTsl . OE-4) 
NOTRAN— logical  VARIABLE  WHICH  IF  .TRUE.  SUPPRESSES 

calling  of  the  TRANSITION  SUBROUTINE  FOR  LAMINAR 
BOUNDARY  LAYERS  (DEFAULTS. FALSE. ) 

HTCORR— LOGICAL  VARIABLE  WHICH  IF  .TRUE.  CAUSES  THE  THERMAL 
ENERGY  INTEGRAL  EQUATION  TO  8F  SOLVED  FOR  THE  HEAT 
FLUX  CORRECTION  PARAMETER,  CHT  f DEFAULTS. FALSE. ) 
ERROR— LOGICAL  VARIABLE  WHICH  IF  .TRUE.  CAUSES  INTERMEDIATE 
H  VALUES  AND  VARIABLES  ASSOCIATED  mITH  THE  TURBULENT 
DISSIPATION  integral  AND  HEAT  TRANSFER  CORRECTION 
parameter  to  be  printed  for  DEBUGGlNf-  PURPOSES 
(DEFAULTS. FALSE.  ) 


Figure  2.  Explanation  of  variables  for  question  12. 
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Sample  Card  Deck 


N2/O2  Case 


0.967  20.0 

2 

N2  0.79 

02  0.21 

300.00 

GSRS  BOUNDARY  LAYER 
TURBULENT  I AXISYM  I  MACH  I  THETA 


$BL  HSTART  =1.6$ 

0.0  1 1.888  1.0 


1.888 


Whether  run  on  air  or  not 
2  -  Not  air 

Mach  Number  -  Pressure 
Number  of  elements  in  mixture 
Element/Mole  Fraction 

Temperature  cf  mixture 
Title  of  Run 

Flow,  geometry,  velocity  input  data,  wall 
temperature  data 
Namelist  BL 
Local  data  cards 


6.63  4.1  2.696  4.1 


XI,  YI,  MDI,  RI,  THETI 
NOTE; 

(Same  as  Air  Case  above) 


1  5 


11  21 


31  41 


XI— Axial  Location,  [L] 

YI - Normal  Location,  [L] 

MDI - Freestream  Mach  Number  (or  MSTAR) 

RI - Cross-Sectional  Radius  of  Axisymmetric  Bodies,  or 

Location  Normal  to  Centerline  for  2-D  Bodies  [L] 

THETI - Wall  Temperature  Ratio  (Tadwall-Twall)/(Tadwall-Tstream)  or  Twall/TINF  - 

VII.  SAMPLE  OUTPUT 

A.  Sample  Case  Using  Air 

B.  Sample  Case  Using 
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'tv 

«V 

'.V 


w:--. 


ii 


|4V 


^  2 


cn 


^  0-)  ^  CM  CO  iH 


I  1 
M  U 
O  O 
<f  <r 
00  o 

ON  ON 
m  ON 
CM  ON 

o  o 

H  CM 

«  • 

o  o 

I 


I  I 
W  M 
o  o 

lO 
l/i  CM 
CO 
CM 

CM  m 
m  CM 
so  ^ 


I  I 
Cz3  (U 
O  O 
sO 
<r  so 

r-l  lO 

o  o 

ON  c^ 
CM  CM 


o  o 
I  t 


o  o 

iH  iH 

1  S 
W  M 

o  o 

CM  O 

m  so 
so  m 

Os 

Mf 

00  so 

rM  iTi 


O  ON 

O 


I 

M  pzq 
o  o 
o  r^ 
Mf  tr\ 

O  iH 
00  CM 

a\  to 

ON  so 


o  o  o  o 
I  I 


Os 

o  o 

I  I 

W 

o  o 

rM  O' 
00  iH 

-o 
m  O' 
CO  O' 
30  r-M 


o  o 


o  o 

I  I 
u  » 
o  o 

f''  <J\ 

tn  eg 
Cl  m 
m  o 
vo  eg 
sr 
r~  m 
•  * 
o  o 

I 


'O  in 

o  o 

I  I 
U  M 
O  O 
00  O 
eg  r-( 

n  o 
ei  sr 
eg  eg 
m 
in  eg 
•  « 
o  o 


vO  sr 

o  o 

M  U 
O  O 
m  r» 
O  (7t 
iH  in 
o  tn 
en  00 

vO  VO 

ov  eg 

•  • 

o  o 

I  I 


u 

fS 

<3 

H 


c^ 

O  O 

I  I 

m  M 
o  o 

VO  00 
T-f 
iH  o 
VO  o 

so  00 

O  M 
rH  CM 


CM  CM 

o  o 

I  I 
U  M 
O  O 

o 

00  in 

iH 

m  00 
*H  o 
m  CM 


CM  iH 

o  o 

I  I 

U  M 
o  o 

O  CO 
O'  CM 

C3N  ON 
00 
CO  CM 
CM 

CM  CM 


o  o  o  o  o  o 


?? 

Ui  u 

o  o 

00  CM 

vO 

CM  O 
O  CO 

in  o 
<f  in 
CM  CM 


?? 

Ui  u 

o  o 

O'  VO 
rH  CM 
CO  00 
vD  •O 
O' 

00  VO 
CM  CO 


o  o 
+  + 
U  Ci3 
O  O 

00  f-! 

«H  O' 
O'  VO 
iH  r>. 
00  vD 
O'  CM 
f— i 


o  o  o  o  o  o 


^  ft 


ARRAY  COEF.  ELEMENTAL  CONSTANT  PRESSURE  SPECIFIC  HEAT  COEFFICIENTS 


s  Si  g  S;  2  «  co‘  2  2 


o^vOfOMOoLTot^  C-I.r,cn 

SI§£S2gs|g§s 


n'o^iS 


o^^c^NvorHcooco^om 
^ofCS!5  5Sr^[?;  oooot'^ 

iH 

sSS§2Sig22i  g'^S 

^  vr  VO  00  g  g  gjsg 

iH 

Oooiovor^c^o^^  i^Or^ 

5^  S  5  g  g  S  S  [S  -S 

>-^<rvD,.25^SSSiS;^:g 


r^<r^iAcNoo^^^  ocvcm 

o  g  ^  Si  ;:^  2  ® 

^srvf.r.  cs5SS2 


'oovovcor^ooocor-  moc 

S  5  2  S  f:  5  9  S  ^ 
^•^voSS^SsS 


^■^cscocor^Ocoo  u-iovn 
c]S222fvi'^‘^'^  ovu-lsr 

'H  ^  .n  S  2  5!  S  ^  S!  g  S!  g 


^(NivCoOOvinvDcMcsl  00<to 
00^00vO00CMvOo2  2rvj(v! 
'-iro,Tv3-,-Hcs)r<i»3-r-<  rvif>~,  ^ 


cnovvo 

'3'0O<f 


'0<Nr-.or-vocinoi  lomcvj 

SiHSsssftgss 


>n^<N^^oooinr«  stovo 

‘aS§S2S9^=^  '3-(^Ov 
hcsm2Ih2  2S^£^  fOHoo 

^‘'^'^'NCOcnrH  cvcoro 


■<f«;rvO{ncrir^(y,^QQ 

5  S  2  S  3  S  «  S 
^<n25SsIs^2S  S!22 

)S  S:  'z!  o  2  rH  2  vr 

■-ICNn^  (SI  cn  CM  (NO 


'Ct'CrrOou^^^OrM  cjvoin 

o  m  2  2?  2  2  2 

S222'"3S2'“  oovvo 

VNIIVJ-  CMCMC»1 


vOvf-crincvro^vDo  cmcmio 

"SS|2g§-'-‘-‘S2 

c'jc')<r  ,_icMn  ,_jpgp^ 


lOc^^rccMr'.vocMr^  Of-co 

222SS:'^'^‘^2  2  2o 

225^'^'^'^*^  '■'i  r^r^in 

c^  m  n  r-i  cvj  pT)  p,,  _ 


ELEMENTAL  THERMAL  CONDUCTIVITIES 


VIII.  RESULTS  AND  CONCLUSIONS 


A  boundary  layer  program  has  been  modified  to  account  for  arbitrary 
chemical  composition  flow  streams.  This  program  can  be  run  in  two  different 
modes: 


1,,  Interactively,  on  the  INTERDi  TA  8/32  for  design  calculations. 

2.  Batch,  on  the  CYBER  for  multiple  runs  and  in  conjunction  with 
other  codes. 

A  sample  case  is  given  in  the  Appendix  which  is  the  same  case  as  given 
in  the  original  COMPEL  documentation.  The  difference  is  that  this  analysis 
specifies  the  chemical  composition  of  air  in  orxler  to  indicate  the  new 
capabilities  added  to  the  code. 

These  improvements  add  the  capability  of  analyzing  internal  flows  for 
aeropropuls ion  applications  such  as  nozzles  and  diffusers. 
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appendix  a 

LINEAR  LEAST  SQUARES  FIT 
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For  a  given  linear  function  y(x)  given  by  (Al) 

y  =  a  +  bx  (Al) 

The  residuals  are  given  by  V]^ 
vj^  =  a  +  bxi  -  yi 

V£  =  a  +  bX2  -  yz  (A2) 


Vn  =  a  +  bxn  -  yn 

In  order  to  determine  the  best  fit,  the  sum  of  the  squares  of  the  residuals 
is  a  minimum 

E(a+bx-y)2  =  (a+bxj^-yj^)^  +  (a+bx2-y2)^  •  •  •  (a+l'^n  -yn  ) 

=  f  (a,b)  (A^) 

For  f  (a,b)  to  a  minimum 

2(a+bxj^-yi)(l)  +  2(a+bx2-y2)(l)  +..«+  2(a+bxn-yn) (1)  =  0  (A5) 


•^  =  2(a+bxi-yi)x-|^  +  2(a+bx2-y2)x2  +...+  2(a+bXjj-yn)Xj^  =  0  (A6) 

Dividing  (A5)  and  (A6)  by  2 


If  =  (a+bxi-yi)  + 

(a-fbx2-y2)  +..•+  (a+bx^-y^)  =  0 

(A7) 

If  =  (a+bxj^-yi)xi 

4-  (a+bx2-y2)x2  +...+  (a+bx^-y^)Xn  =-  0 

(A8) 

Rewriting  (A7)  and  (A8) 

n 

n 

na  +  £  Xj^b  = 

L  yi 

(A9) 

i=l 

i=l 

n  n 

o  « 

£  Xj^a  +  b  £ 

Xi  =  £  X^y. 

(AlO) 

i=l  i=l 

i=l 

Solving  for  a  and  b  from  (A9)  and  (AlO) 
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APPENDIX  B 


SAMPLE  CASE  -  AIR  VERSUS  N2/O2 


In  ordsr  to  test  the  reliability  of  the  boundary  layer  computations  on  a 
given  gas  mixture,  it  was  necessary  to  run  a  sample  case.  In  Case  I,  the 
boundary  layer  program  was  run  in  its  original  form  (that  is,  based  on  air) 
and  Subroutine  NOAIR  was  not  used.  In  Case  II,  Subroutine  NOAIR  was  used, 
based  on  a  mixture  of  79%  N2  and  21%  O2,  the  accepted  composition  of  air. 

It  was,  therefore,  necessary  to  arrange  all  input  data  for  both  cases  to 
be  as  near  alike  as  possible.  Subroutine  NOAIR  had  been  found  to  work  on 
the  N2/O2  mixture  at  a  temperature  of  300°.  (That  is,  the  mixture  viscosity, 
the  ratio  of  specific  heats,  viscosity  power  law  exponent,  and  the  Prandtl 
number  all  were  correct.)  Therefore,  it  was  easier  to  change  the  inputs  of 
Program  COMPEL  (without  Subroutine  NOAIR) . 

CALCULATION  OF  SAMPLE  DATA: 

Three  variables  in  Program  COMPEL  had  to  be  calculated  so  that  the  program 
could  be  tested.  These  were  MINF,  RETNFL,  and  G. 

To  get  MINF,  the  equation  MINF  =  u/a  was  used.  An  arbitrary  u  was  chosen 
to  be  1100  ft/sec  (335.28  m/s).  The  value  for  a  was  calculated  following 

a 
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